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Abstract 

We consider the running of the neutrino mass matrix in the Standard Model and 
the Minimal Supersymmetric Standard Model, extended by heavy singlet Majorana 
neutrinos. Unlike previous studies, we do not assume that all of the heavy mass 
eigenvalues are degenerate. This leads to various effective theories when the heavy 
degrees of freedom are integrated out successively. We calculate the Renormalization 
Group Equations that govern the evolution of the neutrino mass matrix in these 
effective theories. We show that an appropriate treatment of the singlet mass scales 
can yield a substantially different result compared to integrating out the singlets at 
a common intermediate scale. 
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1 Introduction 



The discovery of neutrino masses requires an extension of the Standard Model 
(SM) or the Minimal Supersymmetric Standard Model (MSSM), which may 
involve right-handed neutrinos, or more generally gauge singlets. Since there 
are no protective symmetries, these singlets are usually expected to have huge 
explicit (Majorana) masses. This leads to the see-saw mechanism [1], which 
provides a convincing explanation for small neutrino masses. This scenario can 
be realized in many Grand Unified Theories (GUT's) and their supersymmet- 
ric counterparts. For instance, left-right symmetric models and SO(10) GUT's 
include singlet neutrinos, which can get huge masses in several ways, e.g. by 
a Higgs in a suitable representation or radiatively. Furthermore, additional 
singlets may exist, which can also be involved in the see-saw mechanism. 

It is often assumed that all heavy singlet mass eigenvalues are degenerate. 
However, in all the models a large hierarchy of the singlet masses is possible. 
Note that such a hierarchical spectrum may even show up if all elements of 
the singlet mass matrix are of the same order. Democratic mass matrices, 
where this is the case due to discrete symmetries, are an example. Another 
argument for a non-degenerate spectrum follows from assuming a neutrino 
Yukawa matrix Y v which is proportional to the diagonalized charged lepton 
Yukawa matrix Y e , i.e. the relation Y u = c u Y e pa c v diag(lCT 2 , 1CT 3 , 1CT 5 ) holds 
with a constant real number c v . If the neutrino masses are degenerate and of 
the order 1 eV, the see-saw relation k = -4—M u = 2Yj r M~ 1 Y u for the neutrino 

"ew 

mass matrix M v allows to determine the singlet mass matrix M. Mixings 
do not significantly alter this picture, since e.g. bimaximal mixing can be 
accomplished by small modifications of a degenerate M v of the order 10 -2 eV 
or 1CT 3 eV, respectively. Taking for example c v = 100, the mass eigenvalues 
of M are of the order 10 7 GeV, 10 11 GeV and 10 13 GeV for the case at hand. 
It is therefore conceivable that there may be an even larger hierarchy in M 
than in the charged lepton Yukawa matrices. Altogether, there are thus good 
reasons to study the effects of a non-degenerate or even hierarchical singlet 
mass spectrum. 

In this paper, we calculate the Renormalization Group Equations (RGE's) 
for the evolution of the neutrino mass matrix from the GUT scale to the 
electroweak or SUSY breaking scale. We consider the case where the SM and 
the MSSM are extended by an arbitrary number of heavy singlets which have 
explicit (Majorana) masses with a non-degenerate spectrum. Hence, to study 
the RG evolution of neutrino masses several Effective Field Theories (EFT's) , 
with the singlets partly integrated out, have to be taken into account. Below 
the lowest mass threshold, the neutrino mass matrix is given by the effective 
dimension 5 neutrino mass operator in the SM or MSSM, respectively. The 
corresponding RGE's were derived in [2-6]. 
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2 Effective Theories from Integrating Out Singlet Neutrinos 

Consider the SM or the MSSM with uq additional sterile neutrinos. The 
eigenvalues of the mass matrix M, i.e. the masses of the mass eigenstates 
{iV R , . . . , N^ G }, have a certain spectrum, M x < M 2 < ■ ■ ■ < M na . We will 
consider the general case that this spectrum is non-degenerate. Successively 
integrating out the heavy sterile neutrinos at the thresholds Mj results in 
effective theories, valid in certain energy ranges as depicted in figure 1. 




Fig. 1. Illustration of the ranges of the different theories. The EFT's emerge from 
successively integrating out the heavy fields. "EFT 1" corresponds to the SM or 
MSSM with additional dimension 5 mass operators for neutrino masses. "Full The- 
ory" refers to the SM or MSSM, extended by no gauge singlets. The meaning of 
the variables is explained in the text. 

Before we calculate the RGE's in the various theories, let us specify the mod- 
ifications in the Lagrangians due to the appearance of the heavy neutrinos. In 
the SM above the highest mass threshold ("Full Theory" in figure 1), the ki- 
netic and mass term as well as the Yukawa interaction for the singlet neutrinos 
iV R , i G {1, . . . , nc}, are added: 

^ N = N&vfdjNi + {—NkM^ - (Xu)iMM + h - c ) . (!) 

where iV R := (iV R ) c is the charge conjugate of iV R . / 6 {1, . . . , rip} are flavour 
indices, are the SU(2) L -doublets of leptons, (/) is the Higgs doublet, and 
:= ia 2 (p*. Summation over repeated indices is implied throughout the paper. 
For the calculation of the RGE's, we will work in a basis in which the Majorana 
mass matrix M is diagonal. 
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In the MSSM, the additional gauge singlet Weyl spinors v , which correspond 
to the right-handed Dirac spinors N^, and their superpartners are components 
of the chiral superfields iv c \ The terms of the superpotential containing these 
superfields are 

W m = \^M^ + {Y v ) lf ^\e T n + h.c. , (2) 

where 1^ and [M 2 ) are the chiral superfields that contain the leptonic SU(2) L - 
doublets and the Higgs doublet with weak hypercharge +|. e is the totally 
antisymmetric tensor in 2 dimensions, and a,b,c,d G {1,2} are SU(2) indices. 

The Higgs doublet superfield lb/ 1 ) with weak hypercharge — ^ is involved in 
the Yukawa couplings of the SU(2) L -singlet superfields e c and d c containing 
the charged leptons and down-type quarks, whereas couples to v c and the 
superfield tu c containing the up-type quarks. The part of the superpotential 
describing the remaining Yukawa interactions is given by 

+ (Y d ) gf ^h^e a %( + (Y u ) gf v c ^\e T r\[ , (3) 
where q is the quark doublet superfield. The field content of the superfields is 



|/ = £f + ^/2 0£f + 00 p{ , (4a) 

e Cg = gC 3 + ^2 e Cg + 00 pg ? ( 4b ) 

= V Cj + V2 6v Cj + 99 Fl , (4c) 

q / = q f + V2 6q f + 99 F q f , (4d) 

y C 9 = ~Cg + ^2 u Cg + ee pg ^ ^ 

d c 9 = jcg + ^ 0d C9 + 99 F 9 d , (4f ) 

= + V2 9ftV + 00 FhW t ( 4g ) 

h (2) = 0(2) + ^00(2) + ^ Ffc(a) . ( 4h ) 



By integrating out all singlet neutrinos of the extended SM, one obtains the 
dimension 5 operator that gives Majorana masses to the light neutrinos, 

^ = \* Bf Wl^fa € b e ba <Pa + h.c. . (5) 
The corresponding expression in the MSSM is the F-term of 

K MSSM = -\* a , ie ba hW + h.c. . (6) 

In the intermediate region between the (n — l)th and the nth threshold, the 
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singlets {N^, . . . , N^°} or singlet superfields {a/ 



Cn 



} are integrated 



out, leading to an effective operator of the type (5) or (6) with coupling con- 
stant K g f, where n g f is identical to K g f. In this region, the Yukawa matrix for 

the remaining singlet neutrinos is a (n—1) x rip matrix and will be referred 

(„) 

to as Y u , 



(n)i,i 

\Xv)n-\,\ 



\Yv)n— l,np 











V 



(n) 

--:Y„ 



) 



(7) 



no — n+ 1 heavy, sterile 
neutrinos integrated out . 



The tree-level matching condition for the effective coupling constant at the 



threshold corresponding to the largest eigenvalue M n of M is given by 



(n) 

K 9f 



M„ 



(l+l) 

K 9f 



Air, 



/(n+l)T\ /(n+1) \ 

+ 2\Y V ) gn M-\Y v ) nSM (no sum over n). (8) 



To determine the RGE's, we first calculate the relevant counterterms for the 
effective theories. We use dimensional regularization (with d := 4 — e dimen- 
sions) and the MS renormalization scheme. The renormalization constants 



(") j n ) 

below the nth threshold are denoted by Z, ok, etc., analogous to our notation 
for the coupling constants. 



3 Calculation of the Counterterms 



(n) (n) 

For the one-loop wavefunction renormalization constants Z := 1 + 5Z between 
the thresholds in the extended SM, we find in gauge for U(1) Y and SU(2) L 



(n) 

5Z eL 

(n) 
SZj, 



(n) 
SZ N 



1 



16vr 2 L 
1 



(n) (n) | Q 

YX + YX+ 2^1 + 



(9a) 



16tt 2 



16tt 2 



/(n) (n) \ 

2 Tr [yX) + 2 Tr(y e ty e ) + 6 Tr (Y^Y U ) + 6 Tr (Y}Y d ) 



^B-3)g 2 1 + ^w 



3)^ 2 2 



(n) (n) . 

2 



(9b) 
(9c) 
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For the vertex renormalization constants we obtain 



(n) 



1 



16tt 2 
1 



' (") 4. 1 n (n) 3 o 



$ = --^ 2 (y e t r e ) T( « ) + 2 ( K ) (F e ty e ) - a ( « } 

l07T z L 



(n) 

5M 



o, 



(10a) 

(10b) 
(10c) 



where A is the scalar quartic coupling appearing in the interaction term 
— |A(<^0) 2 - The above quantities are defined by the counterterms for the 
mass and the Yukawa vertex of the sterile neutrinos as well as the one for 
the effective vertex, 



(«) 1 (n) p . 

^mass(N) = — N^SM^ + h.C. , 

(n) 

C ^Y V — ~ 

(„) 



(^) i/ J\^^/ +h . c . 



€=^/*£e°fyd 4^>a + h.C. , 



(11a) 

(lib) 
(11c) 



where the sums over i and j run from 1 to n — 1 . 



In the extended MSSM, only wavefunction renormalization is required except 
for the contributions from the gauge boson - matter interactions. Fixing the 
i?£ gauges and using Wess Zumino (WZ) gauge breaks supersymmetry explic- 
itly, and thus the non- renormalization theorem is not manifest. Hence, the 
counterterms for the vertices do not vanish in general. We use the same nota- 
tion for them as in the SM. The relevant diagrams for the renormalization of 

(n) 

the the /t-vertex are the gauge contributions similar to those of the SM, the 
gaugino contributions (figure 2 (a)-(d)) and the diagrams from the D-terms 
(figure 2 (e)-(f)). The resulting wavefunction renormalization constants are 
given by 



6Z t 



(n) 

5Z, 



0(2) = 



(n) 

SZ N 



16tt 2 
1 

16^2 



16tt 2 



2 yX + 2 yX + g (& - l )al + 2 



2 Tr (yj-yj + 6 Tr(YXu) 
+ \{t,B + i)gl + \{£,w + ±)gl 



(n)(n), 

4^ 



, (12a) 

(12b) 
(12c) 



and the vertex renormalization constants are 
1 



(n) 



lQn 2 



1 (n) 3 (n)' 

-^ B + 2)g 2 1 Y u + -(Cw + 2)g 2 2 Y u 



(13a) 
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(a) 




^3 A 



(c) 




If' 



(e) 



(2) 



A( 2 ) ^ 



(2) 



it 



(f) 



Fig. 2. Figures (a)-(d) are the contributions from the gauginos A A to the renor- 
malization of the dimension 5 operator in the MSSM. Figures (e) and (f) show the 
D-term contributions. The gray arrow indicates the fermion flow as defined in [7]. 



An) 
OK 



(n) 

SM = . 



167T 2 



(^ + 2)3^+3(^ + 2)9^ 



2 (n) 



(13b) 
(13c) 



4 Beta-Functions in the Effective Theories 



4-1 Standard Model with Additional Majorana Neutrinos 

Using the counterterms calculated in the previous section, we find in the SM 

(») , (n ) 

the following /3-functions (3 K = £tj-/%f for the effective vertex below the nth 
threshold: 

+ 2 Ti(Y}Y e )^+ 2 Tr (y^k + 6 Tt(Y^Y u )k 

+ 6Tr(YjY d fK-3g 2 2 ^+^. (14) 
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The method used to calculate /3-functions from counterterms in MS-like renor- 
malization schemes for tensorial quantities is described in [4]. For the Yukawa 



(n) 



matrix, the /^-function f3y v (n > 1) is given by 



» 



(n) 



3 /("),(«) \ 3 , /(«),(«) \ , 



+ 3 Tr(y u ty u ) + 3 Tr(y d V d ) - -a\ - -a\ 



(15) 



4^ ^ 

Calculating the /3-function for the Majorana mass matrix of the singlets yields 

f(n) (n) \ (n) (n) /(n) (n)\ T 



(») /(„) (n) \ (n) (n) /(n) (n) 

16tt 2 /? m = [yX) M + M (y„YJ 



(16) 



4.2 MSSM with Additional Singlets 



In the MSSM with additional chiral superfields including sterile neutrinos, the 
/^-function for the effective vertex below the nth threshold is given by 



16*$ = (YXr^^Ye) + W M CH) 
2 Tr CH) ( £+ 6 Tr(y u ty u ) ( ^- 2g^ - 6gl 



+ 



2 <") 
K . 



(n) 



For we obtain 

(n) (n) 

16tt 2 /?^ = n 



(n) (n) /(ra) (n) \ 

3 yjn + Y}Y e + Tr (yjyj + 3 Tr(yjy u ) - g\ - 3g 



and the /5-function for the Majorana mass matrix of the singlets is 

(n) /WWA (n) (") An) 

16tt 2 /? m = 2 (y.yjj m + 2 m (y,yt J . 



(17) 



(18) 



(19) 



The /3-functions for the gauge couplings and for the Yukawa couplings of the 
quarks and charged leptons are not listed here. We found them to be the same 

(n) 

as in the extended SM or MSSM 8 , if one substitutes Y v -> Y v . 



4-3 Calculation of the Low-Energy Effective Neutrino Mass Matrix 



From the above /3-functions, the low-energy effective neutrino mass matrix 
can now be calculated as follows: At the GUT scale, we start with the Yukawa 
matrices Y v and the Majorana mass matrix M for the sterile neutrinos. Using 
the relevant RGE's (15), (16) or (18), (19) (with the superscripts (n) omitted) 
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together with those of the gauge and the other Yukawa couplings, we calculate 
the renormalization group running of Y v , M and the remaining parameters of 
the theory. 

At the first mass threshold, i.e. the largest eigenvalue M na of M, we integrate 
out the heaviest sterile neutrino and perform tree-level matching according to 
equation (8). Note that this procedure is only possible in the mass eigenstate 
basis at the threshold, which is different from the original one at the GUT 
scale, since the RG evolution produces non-zero off-diagonal entries in M. 
Therefore, the mass matrix has to be diagonalized by a unitary transformation, 
M -> U T M U, which leads to the redefinitions N R -> U T N R , v c -> U T v c and 
Y v — > U T Y U of the singlet neutrino fields and their Yukawa matrix.^] 

Integrating out the heaviest neutrino state yields an effective theory valid 
at mass scales below M nQ . The effective dimension 5 operator k that gives 

Majorana masses to the left-handed SM neutrinos appears in this effective 

(°c) („ g) C«g) 

theory. Next, Y u , k, M, Y e etc. are evolved down to the next threshold, the 

(»g) 

largest eigenvalue of the remaining mass matrix M. The RGE's that determine 
the running of the dimension 5 effective operator between the thresholds are 
given by equation (14) or (17), respectively. 

Again, changing to the mass eigenstate basis, integrating out the singlet neu- 
trino corresponding to this threshold and performing tree-level matching gives 
another contribution to the effective dimension 5 operator. The quantities in 
this effective theory are now evolved down to the next threshold and so on. 
This procedure finally yields the low-energy effective neutrino mass matrix. 

4-4 Running of the Mixing Angle in an Example with Two Generations 

Numerical results for the RG evolution of the mixing angle 8 in a generic 
example with two generations of lepton doublets and two singlets are shown 
as solid lines in figure 3 for the SM and in figure 4 for the MSSM. Here, 9 is 
defined as the angle that appears in the leptonic mixing matrix V = UJU U , 
where U e diagonalizes Y^Y e and U v diagonalizes the effective mass matrix of 
the active (non-sterile) neutrinos. Below the lowest threshold, the latter is 
proportional to the coupling k. In the energy region where heavy neutrinos 

are present, the effective Majorana mass matrix of the non-sterile neutrinos is 

f„i (")_ (™) , (») 
given by « + 2YjM~ 1 Y u . 

5 One could worry that the running, which spoils the diagonal structure of M, 
might require a constant re-diagonalization while solving the RGE's, since their 
derivations assume a diagonal mass matrix. However, this is not necessary because 
the RGE's are invariant under the transformations that diagonalize M. 
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M EW Mi M 2 Mam V 



Fig. 3. RG evolution of the mixing angle in the extended SM with 2 generations of 
lepton doublets and 2 singlets. We used Mgut = 10 16 GeV and the initial conditions 
Mi(Mgut) = 10 8 GeV, M 2 (M GU t) = 10 12 GeV for the Majorana masses of the 
heavy neutrinos at this scale. Besides, we chose the initial values of the Yukawa 
coupling matrices Y v (Mgxjt) to be real with (untuned) entries between 0.025 and 
1. Further explanations are given in the text. 




M EW Mi M 2 Mgut A« 



Fig. 4. RG evolution of the mixing angle 9 in the extended MSSM with 2 gen- 
erations of lepton doublets, 2 sing lets and (</> (2) > / (<^ (1) > =: tan/3 = 35 as well 
as MgusY ~ Mew for simplicity. (A moderate change of the SUSY breaking scale 
Msusy does not change the qualitative picture.) The other parameters are the same 
as in the SM case (cf. figure 3). 
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The transitions to the various effective theories at the mass thresholds lead 
to pronounced kinks in the evolution. For comparison, the dotted and dashed 
lines in figures 3 and 4 show the results when both heavy neutrinos are inte- 
grated out at the higher or the lower threshold, respectively. Obviously, this 
produces large deviations from the true evolution, and the correct result need 
not even lie between the two extreme cases. Although this is only shown for 
the SM in our example, the same happens in the MSSM , if suitable initial val- 
ues for the Yukawa couplings are chosen. Consequently, the correct running of 
the mixing angle cannot be reproduced by integrating out all heavy neutrinos 
at some intermediate mass scale M int e [M 1 , M 2 ] in general. 



5 Discussion and Conclusions 

We have calculated the RGE's for the evolution of a see-saw neutrino mass ma- 
trix from the GUT scale to the electroweak scale in an extension of the SM and 
the MSSM by an arbitrary number of gauge singlets with Majorana masses. 
These masses need not be degenerate and can even have a large hierarchy, as 
pointed out in the introduction. At each mass threshold, the corresponding 
sterile fermion is integrated out, which leads to an effective intermediate the- 
ory and affects the RG evolution of the neutrino masses, mixing angles and CP 
phases. To obtain the low-energy neutrino mass matrix from the Yukawa and 
Majorana mass matrices given at the GUT scale, the RGE's for the various 
effective theories have to be solved. In a numerical analysis for two flavours 
and two singlets, we have found that the renormalization group evolution of 
the mixing angle in the case where the heavy degrees of freedom are inte- 
grated out appropriately differs substantially from that in the case where all 
of them are integrated out at a common scale. The correct running can in 
general not even be reproduced by integrating out all heavy neutrinos at some 
intermediate mass scale. Obviously, similar effects exist for the RG evolution 
of all parameters of a given theory, such as mass eigenvalues, mixings and CP 
phases. 
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